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, I discuss elastic dipole-dipole scattering in QCD at high energies, with emphasis on 

1 the relation between Mueller's dipole picture and the Color Glass Condensate, and 

' on the importance of rare fluctuations for the high-energy limit of the S-matrix. 

1 Dipole dipole scattering: Color Dipoles versus Color Glass 

^ I Let me consider the simplest scattering process that one can think of in perturbative 

i QCD — the elastic scattering between two small color dipoles (or 'onia') at zero 

impact parameter, in the center-of-mass (COM) frame and at relative rapidity Y — 
and address a basic question: What is the high energy limit of the ^-matrix for this 
\ collision ? On physical grounds, the answer seems quite clear: Sy — > as F — > cx) 

I (indeed, at sufficiently high energies, the wavefunctions of the two dipoles contain so 

\ many gluons that the probability Sy that no interaction take place in the collision 

r~| ' is very small, 5*^ <C 1, and it vanishes when Y — > oo). But deriving even this simple 

0^1 limit from perturbative QCD turns out to be quite non-trivial, as demonstrated, 

for instance, by the 'small-x problem' of the linear evolution equations: The BFKL 
^ I equation which is supposed to resum the dominant radiative corrections at high 

. energies, predicts a scattering amplitude Ty = 1 — Sy which rises exponentially 

with Y , thus eventually violating the unitarity bound Ty < 1. 

Over the last two decades, several formalisms have been gradually developed 
which go beyond the BFKL evolution by including those effects which restore uni- 
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^ . tarity at high energies. These effects can be generically characterized as 'multiple 

scatterings', but their precise interpretation, and also their mathematical descrip- 
tion, depend crucially upon the choice of a frame: Whereas in the COM frame, 
unitarity corrections start to manifest themselves, when increasing Y, as genuine 
multiple scatterings ("multiple pomeron exchanges", or "pomeron loops") in 
an asymmetric frame in which one of the two dipoles is much faster than the other 
one, and thus carries most of the evolution, unitarity also appears as non-linear 
effects in the gluon distribution of the evolved dipole, or "gluon saturation" |3l4j . 

This can be best appreciated on a simple example: In the single-pomeron- 
exchange, or BFKL, approximation (valid so long as Ty ^ 1) and in the COM 
frame, the scattering amplitude for two identical dipoles of transverse size tq can 
be estimated as [2]: 

T{ro,ro,Y) aln^ro,ro,Y/2), (1) 

where n(ro, ro, Y/2) is the number density of radiated dipoles of size tq in the wave- 
fimction of the parent dipole of size and with rapidity Y/2, and is the scatter- 
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ing amplitude for two elementary dipoles of similar sizes. The BFKL evolution is 
encoded in the dipole number density, which grows roughly like n(ro, rQ,Y) q'^'Y 
with wp = {4:ln2)asNc/iT (the 'BFKL intercept'). Unitarity corrections become im- 
portant when Ty ^ 1, that is, for Y = Yo with 

Yo-—ln^. (2) 

OJp Qfg 

For Y ^ Yq, and in this particular frame, BFKL fails to describe correctly the scat- 
tering amplitude, but it remains a reasonable approximation for the wavefunctions 
of the evolved dipoles. Indeed, the non-linear effects within each wavefunction are 
rather measured by a^n(ro, r-p, F/2) — the scattering amplitude between a given 
dipole and the other dipoles within the same wavefunction — , which is still small, 
of 0{as), when Y ~ Yq. Therefore, saturation effects in the COM frame start to 
manifest themselves only at higher rapidities Y ^ 2Yo. 

But when the same scattering is viewed in the asymmetric frame in which one 
of the dipoles is nearly at rest, then T(ro, ro, Y) ~ n(ro, tq, Y), with n(ro, ro, Y) 
referring to the evolved dipole. Again, T becomes of 0(1) for Y ^ Yq, but when 
this happens, the non-linear effects in the wavefunction of the evolved dipole are 
also of 0(1), and must be treated on the same footing as the multiple scattering. 

The language of the example above is that of the Color Dipole Picture (CDP), a 
formalism originally developed by Mueller [2] which cannot accomodate saturation 
effects — the onium wavefunction is rather described by BFKL (together with a 
large- A^c approximation), as a collection of dipoles which evolves through dipole 
splitting — , but which can describe unitarity corrections in the COM frame, as the 
multiple scattering between several pairs of dipoles from the two incoming 'onia'. 
This formalism has been numerically implemented by Salam It also lies at the 
basis of Kovchegov's derivation of a non-linear evolution equation for the scattering 
between a dipole and a large nucleus [H], to which I shall return later. 

More recently, a different formalism has been developed, the Color Class Con- 
densate (CGC) which is specially tailored to describe saturation in the wave- 
function of an energetic hadron. The CGC is the matter made of small-x gluons 
in the high-density environment characteristic of saturation. It is characterized by 
a saturation plateau at relatively low transverse momenta : the gluon modes with 
k± < Qs{Y) have large occupation numbers, of 0{l/asNc), but which increase 
only slowly when increasing y or decreasing k± . The saturation momentum Qs (Y) 
|3] is an intrinsic scale generated by the non-linear dynamics; it increases rapidly 
with y, so at large Y it provides a hard scale for the running of the coupling. 

The mathematical language of the CGC theory is that of classical statistical 
physics : This is a theory for classical color fields endowed with a (functional) 
probability distribution which evolves with Y according to a (functional) Fokker- 
Planck equation — the JIMWLK equation |7I8I9| — in which Y plays the role of 
'time'. The classical fields are generated by color sources (the parent dipole plus 
radiated gluons) at rapidities larger than the rapidity Y of interest, whose dynamics 
is therefore 'frozen' by Lorentz time dilation. When increasing Y in one more step 
{Y Y + dY with agdY 1), a new layer is added to the classical fields in 
longitudinal direction, corresponding to 'integrating out' the emitted gluons with 
rapidities between Y and Y + dY. The evolution thus generates a random walk 
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in the configuration space of color fields, in which the elementary step consists in 
the emission of a small~x gluon in the background of the color fields created in the 
previous steps. In this description, saturation eS'ects enter through the fact that 
the probability for induced gluon radiation is non-linear in the background field 
and cannot exceed one. The BFKL evolution is recovered in the limit where the 
background field is weak, corresponding to low gluon occupation numbers. 

Computing a scattering cross-section in the CGC formalism is a delicate task: 
Since the color glass is characterized by strong fields, the standard factorization 
schemes for high-energy scattering are not bound to apply, and in general we only 
know how to describe the collision between a CGC and a simple projectile, like 
a (bare) dipole. Specifically, the S'-matrix for elastic CGC-dipole scattering is 
computed in the eikonal approximation as {xj_ and y± are the transverse coordinates 
of the quark and the antiquark which make up the dipole) 



where W{x±) and U{y±) are Wilson lines describing the color precession of the 
quark, or the antiquark, in the color field yl+ = aa of the CGC, and WV[q!] is the 
probability distribution for this field, which obeys the JIMWLK equation alluded to 
above. By using the latter, one can derive an evolution equation for Sy', however, 
this is not a closed equation, but only the first equation in an infinite hierarchy 
originally obtained by Balitsky 10' . Still, a closed equation for Sy can be obtained 
within a mean field approximation (MFA): This is the non-linear equation originally 
derived by Kovchegov .6: (in a different physical context, though, namely for the 
scattering between a dipole and a large nucleus, where the MFA is better under 
control) , and which is generally dubbed as the Balitsky-Kovcliegov (BK) equation. 

Returning to our original problem of the dipole-dipole scattering, one sees that, 
within the CGC formalism, this scattering is most simply described in an asym- 
metric frame, in which one of the dipoles carries most of the total rapidity and has 
evolved into a CGC, while the other dipole is rather slow and can be described as 
a bare qq pair, without additional gluons. This feature complicates the comparison 
with the CDF formalism, where the same problem is most naturally formulated in 
the COM frame. Still, as shown recently ^Jj through analytic manipulations, the 
two formalisms are in fact equivalent, within the range in Y in which they arc both 
supposed to apply. This equivalence is the first point that I would like to slightly 
elaborate on in what follow. The other point is the role of rare fluctuations in the 
approach of the S'-matrix towards the 'black-body' limit S — T?. (See also Ref. 
[T^ for a related analysis.) 

2 COM scattering between two color glasses 

The interesting rapidity range for comparing CDF to CGC is Yb ^ ^ < 210 (cf. 
Eq. (O) : for Y <^ Yq, both formalisms reduce to the BFKL approximation (and 
thus are obviously equivalent), while for Y > 21o, CDF fails to apply because of 
saturation effects in the wavefunctions of the incoming dipoles (in COM frame). 
The first step towards establishing the equivalence consists in factorizing the S- 
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matrix for the elastic scattering between two color glasses in the COM frame. The 
second step consists in showing that the JIMWLK evolution of the wavefunction 
of a color dipole reduces to the corresponding BFKL evolution (as implemented in 
CDP) in the weak field approximation and for large Nc- 

Concerning the first step, the factorization proposed in reads as follows: 

Sy = J D[afl] Wy/2[c^r] J D[aL] WyM^l] e'S d^-^^'^i^^W^'^ui^^) ^ (4) 

where the symbols L and R stand for the left-mover and the right-mover, re- 
spectively. The exponential factor is recognized as the coupling between the color 
charge density in one system (e.g., = — V^a^) the color field in the other 
system. It describes multiple (eikonal) scattering in the approximation that the 
individual color sources within each system undergo at most single scattering. (But 
global multiple scattering is still allowed, as the simultaneous scattering of several 
constituents from the two systems.) Eq. |0J is correct for Yq < Y < 2Yq since in 
that range none of the two incoming color glasses is at saturation, and multiple 
scattering of a single giuon is indeed negligible. However, global multiple scattering 
is important, because each system involves a large number of constituents. 

As for the second step, one needs to show that the CGC and CDP descriptions 
of an evolved dipole (the onium) become equivalent with each other when the 
CGC formalism is simplified by using the weak field (or BFKL) approximation 
together with the large-iVc limit. To that aim, we have shown in [TT] that (i) the 
parent dipole can be represented as a color glass, and (ii) its evolution with Y, as 
described by the correspondingly simplified version of the JIMWLK equation, can 
be reformulated as the evolution of a system of dipoles, in agreement with CDP. 

The manipulations in Ref. imply that Eq. Q can be rewritten in the 

form expected in CDP namely, as the S'-matrix for the scattering between two 
systems of dipoles. Schematically, 



00 „ . f N N' 

Sy^ drNPN{Y/2) / drN'PN'iY/2) expi ^oi^j) \ , (5) 



where Pn{Y/2) is the probability density for producing a given configuration of 
N dipoles after a rapidity evolution Y/2 (this depends upon the transverse coor- 
dinates of the dipoles and evolves through dipole splitting according to the BFKL 
kernel |2lllj V the integral dFjv runs over the dipole coordinates, and To{i\j) is the 
elementary scattering amplitude (via two gluon exchange) between the dipole i in 
the first onium and the dipole j in the second one. 

Note that, although the exponential in Eq. (0) looks formally like a Glauber 
approximation (the multiple scattering series is resummed as the exponential of 
minus the amplitude for a single scattering), this exponentiation holds only config- 
uration by configuration. After averaging over all such configurations, the resulting 
5'-matrix differs significantly from the simple exponential of the one pomeron ex- 
change". This difference is particularly pronounced in the high energy regime at 



''The one— pomeron-exchange ampolitude in Eq. is recovered from Eq. ^ the linear term 
in the expansion of the exponential there. 
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y <; Iq, where Sy is very smalh The naive exponentiation of the one pomeron 
exchange in Eq. Q would predict {kq is an unknown factor): 

S-y ~ exp| -Koa2n2(ro,ro,r/2)} with n^{Y/2) e'^'^ , (6) 

whereas the Monte-Carlo calculation of Eq. jSJ by Salam rather yields [^j : 

Sy ~ e"''"'^' with a, = a.Nj-rr, k « 0.72 , (7) 

which, although small, is considerably larger than the naive estimate Under- 
standing this difference brings me to my next point, namely: 

3 On the importance of rare fluctuations at high energies 

Why is the 5-matrix in Eq. ITJ approaching the black-disk limit so slowly ? Af- 
ter all, Eq. lO describes the scattering between two typical configurations in the 
wavefunctions of the incoming onia, by which I mean configurations which involve 
a number of dipoles N close to the average value n{Y/2) ~ Q^f^/"^ and for which 
P/v = 0{l). The huge difference between the estimates JT)) and © suggests that, 
at high energy, Sy is rather dominated by rare configurations, which involve only 
few gluons {N <^ n{Y/2)), and therefore have a very low probability to occur 
(P/v(F/2) ^ 1), but which give a much larger contribution to S simply because 
systems with fewer dipoles have a smaller probability to interact, and thus a larger 
probability 5^ to survive without interactions. (Note that the various configurations 
contribute additively to the 5'-matrix in Eq. ((SJ , so the sum there is dominated by 
those configurations which maximize the product Pn{Y/2)Pn'{Y/2)SnxN' ■) 

At this point, it is interesting to note that a result very similar to that in 
Eq. is obtained from the high-energy limit of the BK equation and also 
from approximate solutions to JIMWLK equation valid deeply at saturation [T3). 
Specifically, BK equation yields |14I12| 

Syiro) ^ e-i'^'^^^-^of , (8) 

where Yq is such that QsiYo) ~ I/^Oj £^nd c « 4.88 is the exponent giving the 
energy dependence of the saturation momentum |3]: Q^{Y) ~ (5^(0)e'^"=^. Now, 
in the asymmetric frame in which one dipole is bare and the other one is highly 
evolved, the BK equation describes the scattering between the bare dipole and the 
typical configuration in the evolved one, which is a CGC with saturation momentum 
QsiY) |15I12| . Then, the discrepancy between Eqs. ((HJ and 1© shows that typical 
configurations play very different roles in different frames, thus illustrating the 
strong sensitivity of the physical picture of the high-energy scattering upon the 
choice of a frame. In particular, the configurations retained by the BK equation in 
the COM frame must be some rare configurations, with only few gluons. But then 
it is legitimate to rise doubts about Eq. (jSjl too : Recall indeed that BK equation is 
obtained after a MFA, which should work reasonably well for typical configurations, 
but not also for the rare ones ! But if rare configurations play such an important 
role at high energies in the COM frame, there is no reason why they should be less 
important in the asymmetric frame. In other terms, one cannot trust any result 
like © or (jSJ, which is obtained by including the typical configurations alone. 
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Unfortunately, there seems to be no systematic way to identify the rare con- 
figurations which dominate the ^-matrix at high energies. As a rough criterion, 
the relevant configurations must involve the maximal number of gluons which can 
still give a contribution to S of order one: indeed, further increasing the number 
of gluons would rapidly decrease S, whereas reducing the number of gluons even 
stronger would suppress the probability of the configuration without significantly 
enhancing its contribution to 5". With this criterion in mind, we have been able 
to 'guess' some optimal configurations and then check that they do indeed 
a better job than the configurations retained by the BK equation, in the sense of 
giving a larger contribution to S. Specifically, our best configurations yield 

5^(ro) ^ e-f«=(^-^°)\ (9) 

where as compared to Eq. Q the exponent is now reduced by a factor of two. (This 
result has been confirmed in ^^■) The optimal configurations look differently in 
different frames, but they are rare in any frame. In particular, in the asymmetric 
frame where typical configurations would lead to the result the optimal, rare, 
configurations are those in which the fast onium has evolved into a CGC having a 
lower than normal saturation momentum: Qs{{Y + Yq)/2) instead of Qs{Y)- 

To conclude, let me notice that the exponent in Eq. (O is still larger than the 
one reported in the numerical calculation in Ref. |S] (cf. Eq. ||7J)): c/4 « 1.22 rather 
than K ~ 0.72. So, it would be worth redoing the numerical analysis in order to 
understand the origin of this discrepancy. 
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